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The  current  associated  with  field  emission  is  greatly  dependent  on  the  electric  field  at  the  emitting 
electrode.  This  field  is  a  combination  of  the  electric  field  in  vacuum  and  the  space  charge  created  by 
the  current.  The  latter  becomes  more  important  as  the  current  density  increases.  Here,  a  study  is 
performed  using  a  modified  classical  one  dimensional  (ID)  Child-Langmuir  description  that  allows 
for  exact  solutions  in  order  to  characterize  the  contributions  due  to  space  charge.  Methods  to 
connect  the  ID  approach  to  an  array  of  periodic  three  dimensional  structures  are  considered. 
©  2010  American  Institute  of  Physics,  [doi:  10. 1063/1.3272690] 


INTRODUCTION 

Space  charge,  or  the  effects  of  Coulomb  interactions  be¬ 
tween  emitted  electrons,  and  its  impact  on  the  generation  and 
transport  of  electron  beams  are  important  in  many  practical 
applications,  particularly  when  beam  brightness  and  modula¬ 
tion  characteristics  are  critical  to  the  operation  of  the  device. 
Consequently,  both  space  charge  and  emittance  dominated 
beams  have  received  intense  study,  particularly  as  they  relate 
to  thermal  and  photoemission  sources  for  amplifiers  and  par¬ 
ticle  accelerators.  In  contrast,  field  emission  relies  on 
quantum  mechanical  tunneling  through  barriers  thinned 
through  the  application  of  fields  on  the  order  of  10  GV/m  at 
the  emission  site.  As  the  highest  macroscopic  fields  that  can 
be  generated  are  smaller  (100  MV/m  for  rf  photoinjectors, 
1-50  MV/m  for  vacuum  electronic  devices,  etc.)  field  en¬ 
hancement  must  be  utilized  by  making  the  emission  sites 
have  sharp  features.  ’  While  field  emitters  that  resemble 
pyramids,  ellipsoids,  or  whiskers®  in  microfabricated  arrays'® 
boast  current  densities  at  the  emission  sites  many  orders  of 
magnitude  larger  than  for  thermal  and  photoemission 
sources,  the  emission  area  is  generally  on  the  order  of 
100  nm^  per  site.  The  high  curvature  of  the  site  geometry 
makes  modeling  of  field  emitters  in  modern  particle-in-cell 
(PIC)  codes  an  issue:  simple  models  of  space  charge  or  emit¬ 
tance  from  these  sources  are  unable  to  fully  exploit  the  one 
dimensional  (ID)  techniques  often  used  in  the  analysis  of 
other  electron  sources  or  treat  many  separate  emitters  acting 
in  concert.  Nevertheless,  there  have  been  several  theoretical 
efforts:"  in  particular,  the  treatment  of  the  ID  field 
emission/space  charge  problem  (Ref.  17  and  references 
therein),  the  usage  of  PIC  codes  to  treat  field  emission  (Ref. 
18  and  references  therein)  and  the  discussion  of  numerical 
issues  involved  in  modeling  space-charge-limited  flow  (even 
in  the  ID  case  using  PIC  codes  can  lead  to  unusual  behavior 
as  a  consequence  of  cell-size  related  problems),'®  that  are 
pertinent. 

The  present  study  is  the  first  of  a  two  part  effort  to  de¬ 
scribe  the  coupling  of  ID  and  three  dimensional  (3D)  ap¬ 


proaches  to  field  emission  in  a  manner  that  allows  for  the 
effects  of  space  charge  on  the  operation  of  microfabricated 
structures  to  be  addressed.  Here,  a  methodology  for  the 
analysis  of  space  charge  forces  on  the  emission  process  is 
described,  and  a  manner  in  which  ID  methods  may  be 
brought  to  bear  on  arrays  of  emitters  operating  together  is 
presented  (the  second,  and  separate,  study  concerns  the  im¬ 
pact  of  space  charge  on  individual  3D  structures  and  builds 
on  methods  introduced  in  the  analysis  of  dark  current  and 
emittance.  The  objective  is  to  provide  a  framework  to  in¬ 
vestigate  field  emitters  without  intensive  numerical  efforts  in 
a  manner  amenable  to  PIC  codes  for  when  space  charge  is 
an  issue  and  the  cathode  area  is  but  a  small  region  of  the 
simulation.  Applications  that  rely  on  field  emission  will  ben¬ 
efit,  and  such  applications  include  (but  are  not  limited  to) 
electron  beam  lithography  ’  and  transmission  electron 
microscopes,^"'  spacecraft  propulsion,^^’^®  millimeter-wave 
vacuum  electronic  amplifiers  and  terahertz  devices,  ’  and 
particle  accelerators  and  free  electron  lasers  (FELs).^®^^' 


THE  1D  MODEL 


Poisson’s  equation  in  one  dimension  can  be  written 


(1) 


where  q  is  the  elementary  charge,  p  is  a  number  density,  and 
y(x)  is  potential  energy.  The  current  density  J  is  given  by 

J=qvp=qp\ —  ,  (2) 

\  m  / 


where  (as  in  Table  I)  v  is  the  electron  velocity  and  m  is  the 
electron  mass.  The  second  equality  comes  from  the  assump¬ 
tion  that  the  velocity  is  zero  at  the  cathode  where  y=0.  In¬ 
troduce  the  dimensionless  quantities  tp,  y,  and  j  such  that 

V=<pV,, 
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TABLE  I.  Symbols  and  parameters. 


Symbol 

Definition 

Value 

Unit 

Fundamental  constants  and  parameters 

m 

Electron  mass 

510  999 

eV/c2 

c 

Speed  of  light 

299.792 

nm/fs 

h 

Planck’s  constant 

0.658  212 

eV  fs 

4 

Unit  charge 

1 

4 

Rco 

Rydberg  constant 

13.605  7 

eV 

Fine  structure  constant 

1/137.036 

Bohr  radius  fi!  a^jnc 

0.052  917  7 

nm 

eo 

Permittivity  of  free  space 

5.526  35X  10“- 

q^leY  nm 

Q 

a^JidA 

0.359  991 

eV  nm 

F 

Field  at  cathode 

eV/nm 

V 

Potential 

eV 

Copperlike  parameters 

A 

Chemical  potential 

7 

eV 

kp 

13.554  6 

1/nm 

Uf 

Fermi  velocity 

1.569  19 

nm/fs 

Work  function 

4.5 

eV 

FN  field  emission  parameters 

B 

7pN  parameter,  Eq.  (20) 

65.207  3 

eV/nm 

K 

/pN  parameter,  Eq.  (20) 

0.772  808 

v(y) 

Elliptical  integral  function 

l_y2[l_iln(v)] 

>(y) 

Elliptical  integral  function 

1.061 

Vo 

^-1/2 

0.606  531 

x  =  yL, 


(3) 


where  is  the  anode  potential,  is  a  characteristic  current, 
and  L  is  the  distance  between  the  anode  and  the  cathode 
surface.  Equation  (1)  becomes  (compare  to  the  treatments  of 
Refs.  32  and  33) 


dy'^ 


<p{y)  = 


j 


(4) 


where  is  defined  by 

and  7cl  the  familiar  Child-Langmuir  (CL)  maximum  cur¬ 
rent  that  can  be  drawn  across  an  anode-cathode  gap  with  no 
initial  velocity  of  the  electrons  and  zero  field  at  the  cathode. 
Solving  Eq.  (4)  with  appropriate  boundary  conditions  is 
straightforward  and  has  been  done  in  the  literature  (see  Ref. 
17  and  references  therein),  but  the  solution  is  synopsized  in 
the  present  notation  with  an  eye  to  eliminate  problems  with 
the  appearance  of  nonphysical  results.  Representing  the 
product  of  charge  q  and  field  at  the  cathode  by  F=fVJL, 
where  /  is  dimensionless,  solutions  to  Eq.  (4)  are  then  given 
by 


dip 


(6) 


A  second  integration  yields  [compare  Eq.  (3)  of  Ref.  11] 


+f  =  6fy.  (7) 

Invoking  the  boundary  conditions  that  at  y=l,  9=1,  then 
gives 

{Aj+fy'H2j-f)+f  =  6f.  (8) 

Equation  (8)  is  an  exact  equation  for  the  ID  diode  and  does 
not  have  unphysical  solutions:  manipulations  with  it  should 
be  done  with  care  in  order  not  to  introduce  them.  Equation 
(8)  can  be  written  as 

j=l[2  +  i2-3fhfY73fl  (9) 

which  provides  a  universal  relation  between  the  current  and 
the  field  at  the  cathode  independent  of  the  relation  between 
the  field  at  the  cathode  and  the  emitted  current  (see  also  Ref. 
17).  It  encapsulates  the  two  most  familiar  limits:  in  the  limit 
that  the  current  vanishes,  then/=  1,  or  the  field  at  the  cathode 
F=VglL-,  conversely,  if  /=0,  then  y=4/9  (the  CL  current 
density).  Observe  that  the  field  /  at  the  cathode  cannot  ex¬ 
ceed  1  (its  vacuum  value)  when  space  charge  is  absent  and 
there  is  no  screening. 

We  shall  consider  particular  solutions  depending  on  the 
behavior  of  j(f},  namely,  a  linear  and  quadratic  dependence, 
which  can  be  handled  analytically,  and  a  Eowler-Nordheim 
(EN)  linear  field  model. Using  the  image  charge,  or 
Schottky  barrier  lowering,  factor  is  reserved  for  the  3D 
case  and  shall  be  considered  separately.^®  Before  considering 
these  cases,  we  note  that  Eq.  (8)  can  be  manipulated  and 
re-expressed  in  a  convenient  form  without  square  roots  as 

3/(1-/)-7(4-9;)  =  0,  (10) 

which  bears  a  relationship  to  Eq.  (13)  of  Eorbes*^  [note, 
however,  that  Eq.  (10)  allows  unphysical  solutions  as  the 
second  root  of  the  resulting  quadratic  equation,  e.g.,  /=  1  and 
y=4/9,  satisfies  Eq.  (10)  but  is  not  a  solution  of  Eq.  (8)]. 


Linear  current-field  relationship 

Eor  j  =  af,  that  is,  the  current  is  linear  in  the  cathode 
surface  field,  then  according  to  Eq.  (10)  the  field  satisfies 

0  =  3/(l-/)-fl/(4-9fl/).  (11) 

The  physically  relevant  solution  is  given  by 

/=  |{l  H- 3a^H- (1  -  3a)  Vl  +  fa  +  (12) 

which  yields  /=1  when  a  =  0  (and  thus  y  =  0).  The  case  a 
=  1/3  corresponds  to /=2/3  and  j=2l9  (half  of  the  space 
charge  limit).  The  behavior  of  Eq.  (12)  is  shown  in  Eig.  1  in 
a  form  that  shows  the  asymptotic  convergence  of  /  to  the  CL 
limit:  the  y-axis  is  the  ratio  of  j  =  af  with  the  space  charge 
limit  7cl=4/9.  The  increasing  influence  of  space  charge  can 
be  seen  in  the  small  a  limit,  whereas  the  approach  to  space 
charge  limited  current  can  be  seen  in  the  large  a  limit:  these 
limits  are 
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FIG.  1.  Behavior  of  /  as  given  by  Eq.  (12)  for  the  linear  current-field 
relationship  (the  gray  dashed  line  represents  the  CL  limit). 


i(a  1)  = 


a 

1  -h  -h 


j{a  >  1)  « 


(13) 


j{a  <  1)  ^ 


1,5  22’ 

1  -h 


jia  >  1)  =  (16) 


where  Eq.  (16)  suggests  that  at  a=  1,  y —  0.333,  and  0.254  for 
the  small  and  large  expansions,  respectively,  as  compared  to 
the  exact  value  y  =  0.291.  Also,  observe  that  /  is  approxi¬ 
mately  unity  for  small  values  of  a:  specifically,  to  leading 
order  in  a,  then /(a ^  1)=  1 -(4a/3)  for  both  the  linear  and 
quadratic  cases. 

The  question  of  what  values  of  the  parameters  determine 
when  the  current  becomes  space  charge  dominated  is  very 
important  in  practice.  For  the  ID  geometry  and  the  simple 
emission  models  considered,  the  answers  appear  in  Figs.  1 
and  2;  the  current  is  half  of  the  CL  limit  (i.e.,  j  =  2l9)  when 
a=\l3  and  1/2  for  the  linear  and  quadratic  dependences, 
respectively.  These  quantities  can  be  easily  related  to  physi¬ 
cal  units  by  using  Eq.  (3). 


The  expansions  are  reasonably  good:  at  a=l  Eq.  (13)  sug¬ 
gests  that  y  — 0.346  and  0.333  for  the  small  and  large  expan¬ 
sions,  respectively,  as  compared  to  the  exact  value  y  =  0.367, 


Quadratic  current  field 

When  j  =  af^,  then  Eq.  (10)  becomes 

0=/[3(l-/)-fl(4-9fl/)].  (14) 


FN  current-field  relation 


For  the  triangular  barrier,  the  FN  current-density  versus 
field  relationship  Jp^iF)  [i.e.,  Eq.  (8)  of  Ref.  35]  can  be 
expressed  as 


JpuiF) 


(17) 


where 


The  physical  solution  of  Eq.  (14)  is 

/=^[l-(l-2fl)V4fl+l].  (15) 

6a 

The  behavior  of  the  current  j  =  af'  found  from  Eq.  (15)  is 
shown  in  Fig.  2  in  a  form  that  exhibits  its  asymptotic  con¬ 
vergence  to  the  CL  limit  as  done  in  the  linear  case.  As  be¬ 
fore,  the  influence  of  space  charge  in  the  small  and  large  a 
limits  are 


ln(1+a) 


FIG.  2.  Behavior  of  /  as  given  in  Eq.  (15)  for  the  quadratic  current-field 
relationship  (the  gray  dashed  line  represents  the  CL  limit). 


1  /  jURg,  Vq 

7r(yt  -I-  <I>)  V  (J) 


and  terms  are  as  defined  in  Table  I.  The  work  function  is, 
in  the  absence  of  the  image  charge,  equivalent  to  the  height 
of  the  triangular  barrier  above  the  Fermi  energy  /x  for  a 
metal.  The  form  of  Eq.  (18)  makes  the  dimensionless  nature 
of  a  and  b  transparent.  Their  dependence  on  the  anode  po¬ 
tential  Vg  may  give  pause,  but  recall  that  the  current  density 
depends  on  the  held  strength  at  the  cathode  (i.e.,  /)  rather 
than  the  anode  potential,  and  that  y  is  scaled  by  the  CL  result. 
Equation  (8)  can  be  written  as 

6;'2  - 

(19) 

J 

The  denominator  vanishes  for  y=/^/2.  For  the  left  hand  side 
is  finite,  the  numerator  must  vanish,  or/^=(3/2)/*.  Thus,  the 
degenerate  root  in  the  FN  case  is  given  by  /=2/3  and  y 
=  2/9  (half  of  the  CL  limit).  In  mksa  units,  this  will  occur 
when 
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FIG.  3.  Onset  of  space  charge  effects  for  a  high  (4.0  eV)  work  function  for 
the  ID  field  emission  curi'ent  density  assuming  a  1  /im  anode  to  cathode 
sepai'ation,  metal-like  parameters,  and  a  triangular  bamer  [no  Schottky  bar¬ 
rier  lowering — Eq.  (17)].  The  dashed  line  represents  the  FN  equation  with 
F=y^/L;  the  solid  thick  line  therefore  corresponds  to  F=fVolL  and  there¬ 
fore  includes  the  effects  of  space  charge. 


^frans 


(20) 


which  can  be  solved  iteratively.  For  example,  for  L=l  /uum, 
yu,=7  eV,  and  <I>=4  eV,  then  Vtrans  converges  to  21.953  keV 
after  eight  iterations  for  a  starting  guess  of  10  keV.  Fuans 
therefore  serves  as  a  transition  point  between  when  one 
branch  of  the  numerical  solution  of  is  taken  {V 

<Ftrans)  versus  the  other  (F>Vtrans)'  this  affects,  for  ex¬ 
ample,  the  maximum  and  minimum  V  considered  in  a  bisec¬ 
tion  method  to  find  the  current  density,  as  in  the  first  branch 
Ftrans  the  larger  value,  whereas  in  the  second  branch,  it  is 
the  smaller.  For  an  L  of  1  /Ltm  the  impact  of  space  charge  on 
the  current  density  versus  anode  potential  is  shown  in  Fig.  3 
for  metal-like  parameters.  Field  emission  results  are  more 
conventionally  shown  on  a  FN  plot  of  1  /  V  versus 
ln[7pN(^)/^^]>  so  the  data  of  Fig.  3  are  recast  in  that 
form  in  Fig.  4  with  V  in  keV  and  J  in  A/cm^.  The  impact  of 
lowering  the  work  function  to  2.0  and  0.5  eV  is  shown  in 
Figs.  5  and  6,  respectively.  The  departure  from  linearity  on 
such  plots  indicates  the  impact  of  space  charge  forces. 


Current  regimes:  Limiting  cases 
Sma//  current 

Small  values  of  a  in  the  linear,  quadratic,  and  more  gen¬ 
eral  cases  such  as  FN  correspond  to  small  j,  and  therefore. 


FIG.  4.  Same  as  Fig.  3,  but  for  the  data  represented  on  a  traditional  FN  plot. 


FIG.  5.  Same  as  Fig.  4,  but  for  a  middle  (2.0  eV)  work  function.  The  thin 
line  labeled  “CL”  is  the  Child-Langmuir  limit  [see  Eq.  (5)]. 


Eqs.  (13)  and  (16)  for  j  for  more  general  relations  can  be 
obtained  from  an  expansion  of  Eq.  (9)  for  small  (1-/):  to 
second  order 

;(/-!)- 1(1 -/)-^(l-/)2.  (21) 

In  conjunction  with  a  Taylor  expansion  of  j{f)  for  /  near 
unity  given  by 

;(/)=j(i)-(i-/)y'(i),  (22) 


where  prime  indicates  derivative  with  respect  to  argument, 
then  to  leading  order. 


jif)  = 


;•(!) 

1 +  ]/(!) 


;•(!) 

1  +  8  ’ 


(23) 


which  is  equivalent  (to  order  a)  to  Eqs.  (13)  and  (16)  with 
regards  to  the  denominator.  Eor  the  EN  j(f}  then  Eqs.  (17) 
and  (23)  imply 


(24) 


and  arises  from  increasing  space  charge  when  the  cathode 
efficiency  (larger  a  and  smaller  b)  is  larger. 

The  correction  factor  Sp^  that  arises  from  space  charge  is 
not  difficult  to  examine  experimentally  because  strong  field 
emission  is  not  required.  Importantly,  if  the  anode  potential 
Vg  and  the  anode-cathode  separation  L  are  changed  such  that 
their  ratio  F=VglL  remains  constant,  then  by  virtue  of  the 
definition  of  Sp^  and  j,  it  is  seen 


FIG.  6.  Same  as  Fig.  4,  but  for  a  small  (0.5  eV)  work  function  and  small 
(0.5  eV)  chemical  potential. 
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s^=Ub  +  2)  (25) 

and  so  because  Jchi^o^F/Vg)  scales  as  it  is  seen  that 

SpN  scales  as 


Large  current 


When  the  current  density  is  large  and  j(f}  of  a  more 
complex  form  [e.g.,  Eq.  (17)],  then  numerical  methods  are 
required  to  solve  Eq.  (8).  However,  if  emission  is  so  large 
that  /  is  small,  then  from  Eq.  (9)  it  follows  that 


.  _  (2+V3/+1) 


(26) 


where  [from  Eq.  (5)]yo=4/9.  That  is,  when/>  1/4,  then  j  is 
greater  than  90%  of  the  CL  limit  without  much  regard  to  the 
details  of  the  dependence  of  j  on  /. 

Eor  small  /  and  constant  VJL,  then  the  variation  in  j 
with  a  for  the  linear  (/),  quadratic  {q),  and  EN  relationships 
behave  as 


djg=f, 

(27) 


former  and  u(y)  to  the  latter  (and  will  be  treated  in  greater 
detail  in  the  second  study^®):  the  intention  here  is  to  compare 
the  ID  theory  of  Eq.  (10)  to  the  PIC  simulations  of  Eeng  and 
Verboncoeur,  and  so  their  approximations  shall  be  used  tern- 
porarily.  However,  as  discussed  by  Eorbes,  E^ax 
=  Ufjic  corresponds  to  that  field  which  lowers  the  emis¬ 
sion  barrier  to  the  Eermi  level  (i.e.,  ^=0);  thus,  in  the  usage 
of  the  EN  equation  with  image  charge  correction,  fields 
larger  than  E„ax  disallowed  by  the  approximations  from 
whence  the  EN  equation  is  derived.  If  the  largest  field  con¬ 
sidered  is  10  GV/m,  then  the  work  function  must  be  larger 
than  3.8  eV.  As  a  result,  only  the  findings  summarized  in 
Eig.  11  of  Ref.  18  for  which  the  work  function  is  4  eV  and 
the  anode-cathode  separation  is  1  /rm  are  considered. 

As  shown  elsewhere, the  impact  due  to  changes  in  a 
tunneling  barrier  shape  on  the  exponential  term  in  the  EN 
equation  is  captured  in  the  factor  u(y)  (for  the  image  charge 
potential)  or  analogs  to  it.  That  is,  the  image  charge  barrier 
EN  equation  can  be  made  to  resemble  the  triangular  barrier 
EN  equation  if  the  effective  work  function  depends  on 
the  field  F.  Erom  Eq.  (18),  it  is  seen  that  the  exponential  term 
of  the  triangular  barrier  EN  equation  scales  as 
Equating  the  exponential  arguments  of  the  triangular  barrier 
EN  equation  with  the  image  charge  form,  it  is  seen 

(I,3[2  =  (I,3/2j^(lvafsftcEj.  (29) 


We  note  that  for  small  /  there  is  a  decrease  in  the  magnitude 
of  these  derivatives  which  manifests  iself  in  a  slower  ap¬ 
proach  to  the  CL  limit  (see  Eigs.  1  and  2). 

Comparison  to  PIC  simulations 

PIC  codes  are  widely  used  methods  to  analyze  the  im¬ 
pact  of  space  charge  forces  on  beam  transport.  Electrons  can 
be  launched  into  the  anode-cathode  gap  region,  and  their 
presence  introduces  fields  that  impede  subsequent  emission, 
especially  if  the  emission  is  dependent  on  the  field  at  the 
cathode  surface.  Although  time  dependent  phenomena,  par¬ 
ticularly  in  the  form  of  oscillations  in  the  current  density,  are 
to  be  expected,  such  variations  damp  out  and  a  steady  state 
equilibrium  can  be  approached.  PIC  simulations  addressing 
field  emission  have  been  performed  by  Peng  and 

1  o 

Verboncoeur;  they  examined  the  approach  of  field  emission 
to  space  charge  limited  flow  in  a  ID  simulation.  Here,  the 
findings  of  Peng  and  Verboncoeur  are  compared  to  predic¬ 
tions  based  on  solutions  to  Eq.  (9). 

There  are  differences  between  the  emission  model  used 
by  Peng  and  Verboncoeur  and  the  triangular  barrier  model  of 
the  original  EN  equation:  an  image  charge  corrected  form  of 
the  EN  equation  is  employed  by  the  PIC  simulations  in 
which  the  barrier  is  lowered  by  a  Schottky  factor  by  A<I> 
=  ^afJicF.  Therefore,  the  height  of  the  barrier  above  the 
Eermi  level  p,  is  given  by 

-  iUfJicF  =  {\ -y)^  (28) 

which  defines  the  term  y.  As  discussed  by  Murphy  and 

o 

Good,  the  primary  impact  of  the  Schottky  factor  is  to  alter 
the  a  and  b  factors  in  Eq.  (17)  by  appending  t{y)~^  to  the 


The  coefficient  of  the  exponential  terms  in  both  the  triangu¬ 
lar  and  image  charge  EN  equations  have  different  depen¬ 
dences  on  the  work  function  parameter.  In  keeping  with  the 
“effective”  work  function  concept,  a  factor  of  C  is  appended 
to  the  triangular  EN  equation.  A  naive  equivalence  between 
the  triangular  and  effective  equations  then  suggests 


1 

4<I)f(y)3 


—  (/^+^>eff), 

P  / 


(30) 


which,  for  copperlike  parameters  and  fields  of  3  GV  /  m  (the 
image  charge  barrier  more  closely  resembles  a  triangular  bar¬ 
rier  for  low  field)  is  C  ~  0.42.  Because  the  equivalence  is  not 
exact  (an  effective  work  function  only  roughly  relates  an 
image  charge  barrier  to  a  triangular  barrier)  it  can  only  be 
concluded  that  C  is  on  the  order  of  1/2. 

Solutions  of  Eq.  (10)  using  Eq.  (17)  for  j(f),  and  Eq. 
(18)  for  a  and  b  for  which  is  replaced  by  as  per  Eq. 
(29)  were  calculated.  The  final  j  was  then  scaled  by  the  fac¬ 
tor  C=l/2.  Observe  that  scaling  by  C  does  not  make  the 
scaled  j  [the  solution  of  Eq.  (9)]  the  same  as  the  space  charge 
limited  EN  current  found  by  Peng  and  Verboncoeur;  the  scal¬ 
ing  is  only  to  show  that  the  ID  triangular  barrier  EN  equation 
used  in  Eq.  (9)  anticipates  the  image  charge  EN  equation  in 
that  their  dependence  on  F  is  qualitatively  comparable.  Us¬ 
ing  the  approximation  Peng  and  Verboncoeur  use  of  UFv(y) 
=  1-y*  the  results  of  the  comparison  are  shown  in  Pig.  7; 
the  combination  of  an  effective  work  function  and  a  scaling 
parameter  anticipates  the  PIC  simulations  remarkably  well, 
demonstrating  that  the  ID  solution  using  the  triangular  bar¬ 
rier  EN  equation  is,  in  fact,  a  good  predictor  (albeit  not  a 
replacement)  of  the  behavior  to  be  found  using  the  image 
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FIG.  7.  Comparison  of  the  space  charge  limited  current  as  a  function  of 
applied  field  (ratio  of  anode  voltage  to  anode-cathode  separation)  data  of 
Fig.  11(a)  of  Ref  18  (“Feng  and  Verboncoeur”)  to  Eq.  (10)  using  the  effec¬ 
tive  work  function  of  Eq.  (29)  without  (“unsealed  JFN”)  and  with 
(“scaled  JFN”)  C=l/2:  see  discussion  following  Eq.  (30). 


as  schematically  illustrated  in  Fig.  8.  How  far  from  the  plane 
of  charges  the  anode  in  the  point  charge  model  must  be  be¬ 
fore  the  discrete  nature  of  the  PCM  is  sufficiently  smoothed 
to  approximate  the  planar  cathode  in  the  ID  approach  is  now 
investigated. 

Consider  a  sheet  of  point  charges  spaced  on  a  square 
grid  for  which  the  tip-to-tip  distance  is  (alternately,  the 
pitch  of  the  array)  and  the  magnitude  of  each  charge  is  pro¬ 
portional  to  X:  by  superposition,  the  sheets  due  to  the  other 
X’s  can  be  considered  separately  and  combined  afterward. 
Express  all  lengths  in  units  of  fltt  ^nd  to  take  the  cathode 
surface  to  be  at  z=0  and  the  anode  to  be  at  z=NI2  and  far 
away  (A^S>1).  For  the  anode  potential  to  be  constant  and 
uniform,  image  charges  of  equal  and  opposite  sign  to  the 
cathode  point  charges  are  placed  at  z=N.  The  potential  ev¬ 
erywhere  between  Q<z<NI2  is  then  given  by 


charge  FN  equation  in  PIC  simulations  to  within  a  scale  fac¬ 
tor  C  for  the  space  charge  limited  current. 

The  triangular  barrier  FN  equation  was  used  so  that 
small  work  function  potentials  could  be  considered  without 
running  afoul  of  the  Schottky  factor  lowering  the  emission 
barrier  to  below  the  Fermi  level.  Not  only  for  small  barriers 
but  also  for  barriers,  in  general,  the  j  of  Eq.  (9)  can  be 
suitably  modified  to  accommodate  situations  where  tempera¬ 
ture  and/or  field  make  the  emitted  current  density  complex 
and  possibly  temperature  dependent  (e.g.,  the  general  ther¬ 
mal  field  equation  given  by  Eq.  (36)  of  Ref.  41). 


CO  00 

Vix,y,z)  =  Vo  +  7 - 2  2  <Pjkix,y,z), 

fc=-oo 


<Pjkix,y,z) 


Rc 


R. 


RaRciRa  +  Rc) 


N{N-2z) 

RcRaiRc  +  Rj 


Rc=[<J-x)^  +  {k-y)^  +  z^V’^, 


Connection  of  the  1 D  model  to  an  array 

The  ID  model  regards  emission  of  uniform  sheets  of 
charge  from  a  surface.  Field  emission  arrays,  in  contrast, 
emit  from  a  lattice  of  emitter  sites.  As  is  well  known  in 
electrostatics  (and  explicitly  utilized  by  the  point  charge 
model^°’^*),  conducting  surfaces  can  be  replaced  by  equipo- 
tential  surfaces  at  the  same  potential,  and  the  converse  holds 
as  well:  the  equipotential  surface  signifying  “anode”  in  the 
point  charge  model  becomes  the  “cathode”  of  the  ID  region. 


Rc=[{J-x)^+ik-y)^+iN-z)^V^\  (31) 

where  the  ordering  of  the  arguments  in  q>  reflect  that  the 
potential  satisfies  Vix+n,y  +  m,z)  =  Vix,y ,z)  for  integers  n 
and  m,  and  therefore  x  and  y  may  be  assumed  less  than  1/2. 
For  computational  purposes,  the  last  form  of  <p,j  in  Eq.  (31) 
is  preferred.  Fikewise,  the  z  component  of  the  field  is  given 
by 


FIG.  8.  (Color  online)  Schemmatic  of  an  array  of  emitters  in  the  ungated 
(left)  and  gated  (right)  configurations.  The  ID  approach  applies  for  large  z 
(top)  and  the  unit  cell  approach  for  small  z  (bottom).  An  enlargement  of  the 
gated  emitter  is  shown  on  the  bottom  right,  and  its  Saturn  model  represen¬ 
tation  on  the  bottom  left,  where  the  point  and  ring  charges  are  in  the  gate 
plane. 


f  •  F(x,y,z)  = 


00  00 


7^2  2 

47r£o^2tty=_oo  k=-oo 


‘Pjkix,y,z), 


(Pjkix,y,z) 


N-z 


R 


3 

a 


(32) 


The  other  components  are  analogous,  but  the  z  component  is 
dominant.  How  large  must  z  be  before  the  point  particle  na¬ 
ture  of  the  charges  is  obscured  sufficiently  to  approximate 
the  ID  framework  is  the  question. 

Consider  a  finite  circular  array  of  point  charges  of  diam¬ 
eter  2M  such  that  Eq.  (31)  is  the  M^oo  limit.  It  is  expedient 
to  numerically  analyze  finite  M  and  generalize  to  the  M 
— >  00  case.  For  z  sufficiently  large  that  the  series  can  be  well 
approximated  by  an  integral,  it  follows 
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FIG.  9.  (Color  online)  The  z  component  of  the  gradient  of  the  potential 
along  the  diagonal  defined  by  x=y  for  a  circular  emitters  arranged  on  a 
square  grid  with  (A^,M)  =  (5, 10)  for  the  left  (a),  and  (10,40)  for  the  right  (b) 
for  various  values  of  z  (red  (dashed) =0.4,  blue  (thin  line) =0.8, 
black  (thick  line)  =  1.6).  The  length  is  scaled  by  half  of  the  diagonal  of  the 
total  square  region  (that  is,  \5m).  For  infinite  uniform  charged  parallel 
plates,  the  field  between  them  would  be  constant  (independent  of  z)  and 
equal  to  F^. 


f  •  F  =  lim 

M— >00 


M 

J_  V 

F„  4-tt,- 


2 


M  ^ 


2 

=  1  - 


i  f 

2J0 


zr 


(N-z)f 


(r2  +  z2)3/2  |-^2^(^_^)2]3/2 

z  (JV-z) 


fdr 


2  +  z^  2  ^M^+iN-zf  ’ 


(33) 


where  Fg=\l  e^a^.  It  is  seen  that  the  second  two  terms  in  the 
third  line  of  Eq.  (33)  are  equivalent  to 

z  {N-z)  1 


2VM^  +  z^  2\m^+{N-z) 


rdr(p'jyz)\j2+k2=p- 

J  M 

(34) 


For  Mi>N>z,  then  to  leading  order  the  z  component  of  field 
is  F=Fg[l-{N/2M)].  Clearly,  as  F  approaches  the 

parallel  plate  solution  Fg.  The  finite  M  form  of  Eq.  (32)  can 
be  evaluated  numerically,  and  how  large  z  should  be  can  be 
inferred  from  when  Eq.  (33)  approaches  its  1-{NI2M)  limit. 
Consider  the  variation  in  F^^l Fg  given  by  Eq.  (33)  along  the 
line  x=y,  for  which  variations  are  at  their  maximum,  for 
representative  values  of  {N,M)  of  (5,10)  and  (10,40)  in  Fig. 
9.  For  z  >  1 ,  the  undulations  in  F^  are  largely  absent.  This 
can  also  be  seen  in  the  behavior  of  the  potential  itself  in  Fig. 
10  for  {N,M)  values  of  (4,8)  at  z=0.4  (A)  and  z=  1.6  (B):  the 
color  scales  differ  as  the  potential  increases  with  z,  but  it  is 
clear  that  the  collection  of  pointlike  potentials  for  z  <  1 
merges  into  a  more  uniform,  or  sheetlike,  behavior  for  z 
>  1 .  The  small  values  of  N  and  M  are  to  facilitate  visual 
inspection,  but  field  emitter  arrays  are  characterized  by  far 
larger  M  than  considered  herein,  for  which  F^  is  more  nearly 
constant  within  the  disk  defining  the  active  array  area.  Con¬ 
sequently,  z  values  of  order  unity  are  the  appropriate  loca¬ 
tions  for  which  to  transition  from  a  unit  cell  representation 
(discussed  separately^^)  to  the  ID  approach,  with  the  anode 


FIG.  10.  (Color  online)  Same  as  Fig.  9,  but  showing  the  potential  itself  over 
the  array  for  (A^,M)  =  (4,8)  and  z=0.4  (a)  and  1.6  (b).  The  solid  color  in  the 
center  of  B  shows  that  as  z  increased  beyond  unity,  the  potential  near  the 
center  is  more  uniform. 


of  former  being  the  cathode  plane  of  the  latter.  Issues  remain 
about  the  optimal  choice  of  the  field  of  that  boundary,  and 
how  to  account  for  its  small  but  nevertheless  present  varia¬ 
tion  with  X  and  y,  but  those  questions  shall  be  taken  up  in  a 
separate  study.  For  present  purposes,  the  ID  /  is  related  to 
the  3D  Fg  y'mf^Fg  and 

A  demonstration  of  how  rapidly  the  potential  (or  electric 
field)  directly  above  a  point  charge  converges  with  the  po¬ 
tential  (or  electric  field)  above  the  midpoint  of  four  point 
charges  provides  the  final  indication  of  how  rapidly  the 
ripples  in  Figs.  9  and  10  decline.  As  seen  in  Eq.  (31)  for  the 
potential  and  Eq.  (32)  for  the  field,  it  will  amount  to  finding 
numerically  how  rapidly  the  potentials  off  axis  and  on  axis 
approach  the  same  value.  From  Eq.  (31),  it  is  seen  that  such 
a  question  can  be  most  easily  answered  by  ignoring  Vq  and 
evaluating 


R<piz) 


^,/2+it2<M2^j/s(0,0,z)  -I- 

/l  1  \ 

V’jk\ 


^j2^k2<M2(Pjki0,0,z] 

- Tfrr 

-,-,Z 

Equation  (35)  will  approach  0  as  z  increases.  It  is  seen  that 
adding  ATTBQaffVoIX  to  both  the  numerator  and  denominator 
in  the  fractional  part  (where  the  M— >oo  limits  give  terms 
leading  to  the  potential  and  field  from  and  Rp,  respec¬ 
tively)  will  only  cause  R^  and  Rp  to  decrease  more  rapidly. 
Therefore,  Eq.  (35)  is  a  good  metric  to  consider  how  the 
ripples  fade.  The  terms  and  evaluated  analogously  to 
Eq.  (34),  appear  in  both  the  numerator  and  denominator,  as 
for  large  M,  whether  x=y  =  0  or  1/2  is  of  small  importance. 
We  find 

Sm=^{[M^+{N-  zf]  -  [M^  +  z^]  , 

-  -I  ^ 

"  2[[M^+iN-  zfy^^  [M^  H- 

The  behavior  of  R^z)  and  Rp{z)  are  shown  in  Fig.  11  as  a 
function  of  z  where  M  is  taken  to  be  20  (for  comparison, 
M=10  is  also  shown,  and  A=  10  for  both):  it  therefore  is  a 


+  S’^ 


-  1. 


(35) 
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z-  D 

FIG.  12.  Magnitude  of  the  ripple  amplitude  functions  and  Rf  as  a  func¬ 
tion  of  distance  z-d,  where  z  is  in  units  of  the  tip-to-tip  spacing  of  the  airay, 
for  the  dipole  model  of  Eq.  (37). 


FIG.  11.  Magnitude  of  the  ripple  amplitude  functions  and  If f  as  a  func¬ 
tion  of  distance  z  from  the  cathode  surface,  where  z  is  in  units  of  the 
tip-to-tip  spacing  of  the  array,  for  the  “monopole”  model  of  Eq.  (35). 


measure  of  the  amplitude  of  the  ripples  of  Fig.  9  and  dem¬ 
onstrates  the  rapidity  with  which  those  ripples  vanish.  The 
field  ripples  vanish  less  rapidly  than  the  potential  ripples 
because  (p'  is  more  sensitive  to  variations  than  <p. 

With  a  minor  modification,  it  is  possible  to  consider  a 
cathode  composed  of  sharper  tips  in  a  background  field  F„. 
Let  the  cathode  plane  be  at  z=0  as  before,  but  now  consider 
two  oppositely  charged  particles  symmetrically  placed  about 
the  z  plane  for  each  emitter  which  we  shall  call  the  “dipole” 
model  (it  anticipates  the  dipole  model  that  shall  be  examined 
at  length  in  the  single-tip  space  charge  investigation  consid¬ 
ered  separately),  in  contrast  to  Eq.  (35)  that  uses  only  one 
charge  to  represent  the  tip  (“monopole”  model).  Therefore, 
the  potential  of  the  dipole  in  a  background  field  is  given  by 

00  00 

ydipiix,y,z)  =  F^z  +  — E  E 


1 


1 

Rj 


R±  =  [{j-xf  +  ik-y)^+iz±  (37) 

In  other  words,  the  anode  at  z=NI2  in  the  monopole  model 
has  been  replaced  by  a  background  field  F^.  As  done  with  Vq 
in  the  monopole  model,  the  background  field  Fg  may  be 
neglected  in  the  dipole  analog  of  Eq.  (35),  as  the  ratios  minus 
1  will  be  smaller  with  Fg  than  without  it,  and  so  the  calcu¬ 
lation  without  Fg  serves  as  an  upper  bound.  It  is  seen  that  the 
dipole  model  formulas  for  R^  and  Rp  ctre  obtained  by  the 
replacements  N^2d  and  z^z+d  in  the  monopole  formulas 
for  R  and  S,  but  the  interpretations  are  different;  d  in  the 
dipole  model  is  smaller  than  the  tip-to-tip  spacing,  whereas 
N  in  the  monopole  model  is  much  greater  than  the  tip-to-tip 
spacing.  Performing  the  analogous  calculation  for  the  dipole 
model,  the  results  are  shown  in  Eig.  12  using  M=20. 


Gated  structures 

A  final  complication  is  the  possible  presence  of  a  gate 
near  the  field  emitter,  which  for  field  emitter  array  (EEAs) 


serves  to  modulate  the  array:  “emission-gated”  current  is  re¬ 
quired  by  a  variety  of  technological  applications  from  dis- 
plays  to  microwave  amplifiers  and  EELs.  The  simplest 
model  of  the  gate  (a  metal  plane  collinear  with  the  emitter 
apexes  from  which  disks  centered  about  the  emitter  apex  are 
excised,  as  shown  in  the  expanded  diagram  in  Eig.  8  of  the 
gated  cell)  is  to  add  to  the  charges  that  represent  the  emitter 
a  charged  ring  to  represent  the  gate  (the  so-called  Saturn 
model"^^^).  The  potential  is  then,  in  spherical  coordinates. 


^>saturn('‘> 


47reo| 


I 

r 


(38) 


where  P,(cos  9)  is  a  Legendre  polynomial,  and  where  the 
magnitudes  of  the  charges  ^tip  zas  comparable. 

Therefore  (p  decreases  faster  than  in  the  ungated  case  when 
the  ring  is  absent,  to  leading  order  as  (^tip“?ring)/G  or,  if 
9'tip=9'ring>  2qfjj,apP  for  small  9.  Therefore,  an  array  of 
Saturn  rings  plus  charges  should  coalesce  more  rapidly  to  the 
ID  representation  than  the  ordered  array  of  bare  point 
charges  considered  previously. 

Eor  both  gated  and  ungated  geometries,  an  approach  that 
may  allow  for  the  estimation  of  the  impact  of  space  charge 
on  field  emission  from  arrays  in  a  manner  amenable  to  PIC 
beam  simulation  codes  is  suggested.  PIC  is  used  to  model  the 
injection  and  acceleration  of  charge  bunches  especially  when 
space  charge  complicates  transport  and  causes  emittance 
growth. Indeed,  the  needs  of  PIC  codes  suggest  using 
the  PCM  approach  to  develop  the  cathode  boundary  of  the 
PIC  simulation.  In  so  framing  the  problem,  the  space  charge 
limits  considered  in  the  ID  section  of  this  work  have  direct 
bearing  to  more  comprehensive  and  time  dependent  PIC 
simulations,  thereby  allowing  field  emission  sources  to  be 
treated  by  methods  that  in  the  past  have  been  profitably  used 
primarily  on  planar  cathode  structures. 


UNIVERSAL  FEATURES  OF  FIELD  EMISSION  IN  1D 

In  spite  of  its  simplicity,  the  ID  CL  model  for  a  cathode 
of  unlimited  emissivity  plays  a  very  important  role  in  the 
design  and  analysis  of  various  devices*  as  well  as  for  devel- 
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FIG.  13.  The  ratio  of  9/74,  where  j  is  the  ratio  of  current  density  with  the 
CL  limit,  as  a  function  of  /  (ratio  of  surface  field  F  to  vacuum  field  V„IL) 
for  various  cunent-field  relationships:  L=linear;  2=quadratic,  and  FN 
=Fowler-Nordheim,  for  the  cases  a=b=\. 

oping  effective  computational  schemes  for  more  realistic  ge¬ 
ometries.  The  situation  is  more  complicated  in  the  case  of 
field  emission  because  the  current  density  depends  on  the 
field  at  the  cathode  surface,  but  surprisingly,  Eq.  (9)  allows 
one  to  find  universal  features  of  emission  which  are  indepen¬ 
dent  of  specific  cathode  properties  [as  also  done  by  Forbes,*^ 
who  found  a  simplified  version  in  his  Eq.  (15),  although  Eq. 
(9)  here  does  not  introduce  extraneous  roots].  As  noticed  by 
Forbes,  and  shown  by  Eq.  (9),  j  is  a  function  of  /  only,  and, 
in  particular,  how  close  the  current  produced  by  field  emis¬ 
sion  approaches  the  CL  limit  is  determined  solely  by  the  field 
at  the  cathode  surface  (in  units  of  the  vacuum  field).  To  show 
this,  a  plot  of  9y/4,  with  j  given  by  Eq.  (9),  is  compared  to 
the  example  lines  ji=f,  jq=f',  and  ypN=/^  6xp(-l //)  is 
shown  in  Fig.  13.  The  intersection  points  of  the  former  with 
the  three  latter  lines  give  the  value  of  /  that  solves  Eq.  (9)  for 
the  particular  current-field  relation  (linear,  quadratic,  or  FN, 
respectively). 

As  demonstrated  in  Fig.  12,  9y74>0.9  as  soon  as  / 
<0.27  (alternately,  F<0.27V„IL),  and  even  when /=  1/2, 
the  current  is  not  far  from  the  CL  limit  (9_/74~ 0.697).  The 
conclusion  that  a  field  with  a  magnitude  of  a  quarter  of  the 
vacuum  field  V„IL  is  sufficient  to  make  emission  almost 
equal  to  the  CL  limit  can  be  observed  in  experiments.  When 
9j/4  approaches  unity  and  /  approaches  zero,  then  the  fol¬ 
lowing  relationship  is  obtained  from  Eq.  (9)  for  F^VJL, 

“  ( ^)  VoycLiVo,L)  -  mi  (39) 

where  we  have  returned  to  dimensioned  units  and  use  has 
been  made  of  Eq.  (5).  Thus,  for  a  given  CL  limit,  F  scales  as 

yll4 

CONCLUSION 

The  relation  of  space  charge  to  field  emission  is  an  im¬ 
portant  problem  because  of  the  strong  variation  in  the  emit¬ 
ted  current  with  the  field  that  exists  at  the  emission  site. 
Consequently,  charge  between  anode  and  the  emitter  (or 
emitter-gate)  boundary  bears  a  complex  relation  to  the  volt¬ 
ages  and  separation  distances  defining  the  diode  region. 
Three  studies  of  space  charge  and  its  impact  on  field  emis¬ 
sion  are  therefore  indicated.  First,  a  ID  analysis  of  the  basic 


Poisson  relation  in  a  diode  region  is  required.  Second,  a  rep¬ 
resentation  of  a  3D  field  emitter  structure  adaptable  to  the  ID 
theory  is  needed.  Third,  and  finally,  a  method  to  approximate 
the  boundary  needed  in  the  ID  problem  by  considering  the 
behavior  of  a  3D  periodic  array  is  needed.  In  the  present 
study,  the  first  and  third  issues  have  been  addressed,  and  the 
second  shall  be  reported  separately.  Although  the  present 
analysis  is  steady  state,  its  results  are  commensurate  with  the 
asymptotic  behavior  of  the  PIC  simulations  of  Feng  and 
Verboncoeur.  The  transition  from  3D  to  ID  for  the  case  in 
which  field  emission  arises  from  a  surface  that  has  irregulari¬ 
ties  or  is  purposely  modified  in  such  a  way  that  the  average 
spacing  between  the  irregularities  (spikes)  is  H,  than  at  a 
height  FI  above  the  surface  of  the  cathode  toward  the  anode, 
the  electric  field  is  indistinguishable  from  the  ID  field,  which 
is  of  consequence  when  cell-size  constraints  that  are  already 
important  in  PIC  codes  (e.g..  Ref.  19)  are  an  issue. 
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